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Abstract 

First-principles full potential calculations based on Zeroth Order Regular Approximation 
(ZORA) relativistic Hamiltonian and Kohn-Sham form of Density Functional Theory (KS DFT) 
in local spin density approximation (LSDA) are reported for group IIIA-VA (InAs, GaAs, InP) 
semiconductors. The effects of relativity are elucidated by performing fully relativistic, scalar 
relativistic, and nonrelativistic calculations. Structural and electronic band structure parameters 
are determined including split-off energies, band gaps, and deformation potentials. The nature 
of chemical bonding at the equilibrium and under hydrostatic strain is investigated using pro- 
jected (PDOS) and overlap population weighted density of states (OPWDOS). ZORA results are 



^— s . compared with Augmented Plane Wave plus Local Orbitals method (APW-l-lo), and experiment. 

*0 I Viability and robustness of the ZORA relativistic Hamiltonian for investigation of electronic and 



structural properties of semiconductors is established. 
PACS numbers: 
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I. INTRODUCTION 



There is a great interest in electronic and structural properties of group IIIA- VA materials 
due to their wide spread applications in semiconductor devices. In particular, InAs/GaAs 
and InAs/InP semiconductor quantum dots (QDs) [1] have shown great promise [21] in quan- 
tum information applications such as is the generation of entangled photon pairs (EPPs) on 
demand [sHSl]. 

Atomistic modeling of semiconductor nanostructures may require input from accurate 
Density Functional |7H9| calculations in cases when experimental data is not available. 
Therefore, it is important to understand which material parameters are well reproduced 
with "standard DFT" and this work is a contribution in this area. Three factors determine 
accuracy of Density Functional calculations 1) Model exchange-correlation functional; 2) 
Representation of single-particle orbitals (atomic orbitals, plane-waves, real space grids) and 
representation of ion-electron interaction [ab initio pseudopotentials, full potential schemes); 
3) Treatment of relativity. The assessment of the accuracy of exchange-correlation function- 
als is beyond the scope of this contribution. The main objective of this work is to perform a 
detailed study of structural and electronic structure properties of InAs, InP, and GaAs semi- 
conductors using highly accurate representation of single-particle orbitals and ion-electron 
interaction and to assess the role of relativity in these calculations. 

First-principles calculations on group IIIA-VA semiconductors based on Kohn-Sham form 
of Density Functional Theory 7H9| have already been performed in the past lCHl5l|. The 
computational approach used in these calculations generally evolved from ab initio pseudopo- 
tential calculations to more elaborate full potential (FP) augmentation schemes 16l-ll9l| such 
as Linearized Augmented Plane Wave (LAPW) and Augmented Plane Wave plus local or- 
bitals (APW+lo) methods. LAPW and APW+lo have now become the methods of choice 
when accuracy considerations have the highest priority. 



There is also a less known full potential method due to te Velde and Baerends 20|, l21| 
which makes use of Bloch basis set made up of numerical and Slater type atomic orbitals 
(NAO/STO basis), a basis which is excellent for the accurate representation of the electron 
density. The approach due to te Velde and Baerends is implemented in BAND program 22] 
and is capable of treating chemical elements throughout periodic table using Zeroth Or- 
der Regular Approximation (ZORA) relativistic Hamiltonian 23142 7l|. ZORA approach is 



designed to capture scalar relativistic effects such as s and p "orbital contraction" (stabiliza- 
tion) and d "orbital expansion" (destabilization) as well as spin-orbital splitting for electrons 
with angular momentum / > 0. While ZORA Hamiltonian is very well established among 



quantum chemists, relativistic ZORA calculations on solids are much less common 28N33|. 
especially, in comparison with a large volume of calculations employing LAPW and APW-l-lo 
methods. Therefore, further assessment of ZORA performance in solids is important. 

In this work, 1 will employ ZORA Hamiltonian to perform a detailed study of structural 
and electronic band structure properties of InAs, InP , and GaAs semiconductors in zinc- 
blende phase. The effects of relativity are elucidated by performing three sets of calculations 
1) nonrelativistic 2) scalar relativistic and 3) relativistic with variational treatment of spin- 
orbital coupling (or fully relativistic). ZORA Hamiltonian is applied to calculate electronic 
band structures, band gaps, and deformation potentials. The results obtained with ZORA 
Hamiltonian are compared to those obtained with APW+lo method, and with experiment. 

I perform a detailed analysis of electronic structure of InAs, InP and GaAs in terms of 
Projected and Overlap Population Weighted Density of States (PDOS/OPWDOS analysis). 
Whereas it is a relatively general practice to report and discuss PDOS, the usage OPWDOS 



is much less common. The OPWDOS analysis, popularized 3J] by the Nobel Prize winner 
Roald Hoffmann, provides a clear pictorial representation of bonding, non-bonding, and 
anti-bonding orbital interactions and is deemed useful. I will demonstrate the usage of 
OPWDOS plots and explain how they add to our understanding of chemical bond in InAs, 
InP, and GaAs. 

Atomic units h = e = rrie = I are used throughout unless otherwise specified. 

II. COMPUTATIONAL APPROACH 

The first-principles KS DFT calculations with ZORA Hamiltonian are carried out 



with BAND program 



20y22 



22| . BAND makes use of periodic boundary conditions (PBC) and 



explicit Bloch basis composed of numerical and Slater type atomic orbitals (NAO/STAO 

nn 

basis). I also perform calculations with APW+lo method |16|, llll as implemented in EXCIT- 
ING program 35|]. Both methods are capable of including spin-orbital coupling variationally 
and I choose to do so. To clarify effects of relativity, I also present results of scalar relativistic 
and nonrelativistic calculations. 



The relativistic effects in BAND are treated in the Zeroth Order Regular Approximation 



(ZORA) approach of van Lenthe and co-workers 23l-l27l|. The details of ZORA implemen- 



tation in BAND are given in Ref. 29|. In ZORA approach, the kinetic energy operator is 



replaced by the ZORA expression T^ora 

where p is the "momentum operator" (p = — iV in the absence of a magnetic field); Vsapa{j) 
- is a sum of atomic potentials (SAPA), an approximation to the total effective potential in 
the ZORA kinetic energy operator; and a = {ax, <Jy, cTz} - is a vector made up of the Pauli 
matrices. Introducing the notation 

^ = l-VsAPA{r)/2c^' ^^^ 

the ZORA kinetic energy becomes 

T^°^^ = a-Py<T-p = PyP + ^a-(Virxp). (3) 

In the last equation, T^^^^ was split into the so-called scalar relativistic Tg^^^ and 
spin-orbital T^q^^ terms, where T^q^^ = 1/2(T ■ (VK x p). The nonrelativistic limit can 
be obtained by setting K ^ 1. I will refer to calculations with T^^^^, Tgj^^^, and K = 1 
as (fully) relativistic (ZORA FREL), scalar relativistic (ZORA SREL), and nonrelativistic 
(NREL), respectively. 

In my BAND calculations, I use basis set of triple zeta quality (TZ2P in BAND's nota- 
tion) taken from the program's database. The core states are obtained from the full potential 
atomistic calculations and are kept frozen during the self consistent field (SCF) procedure. 
The valence states are expanded in terms of the NAO/STAO Bloch basis set functions or- 
thogonalized on the core states (VOC basis). The Hamiltonian matrix elements are evaluated 



using highly accurate numerical integration scheme [36|]. The Brillouin zone integration is 
carried out using accurate quadratic tetrahedron method 37|, |38|| with 65 symmetry unique 
k-points spanning the irreducible Brillouin zone (IBZ). The "default" convergence criteria 
are used to terminate the SCF procedure. 



APW-l-lo calculations are carried out with EXCITING program 351]. The local orbitals 
and linearization energies are taken from the program's database. The "core" states are 
treated fully relativistically and self-consistently in the spherical approximation, whereas the 



"valence" states are treated using the second-variational Hamiltonian. The IBZ is sampled 
using {8x8x8} uniform mesh of k-points. 

The exchange-correlation is treated within the local spin density approximation 



E 



(LSDA) I8,l9|. BAND and APW+lo LSDA calculations are carried out using Vosko-Wilk- 
Nusair |39| and Perdew-Wang [40] parameterization of the correlation energy, respectively. 
The calculations are performed using "primitive" face-centered cubic cell with two atoms 
per cell. The zinc-blende crystal structure was assumed. 

The structural parameters are obtained by varying the lattice constant (from -20 to 20% 
of the equilibrium volume) and fitting the total energies to the Murnaghan equation of 
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The electronic structure is analyzed using projected density of states (PDOS) and overlap 
population weighted density of states (OPWDOS). PDOS is defined for a function x^ or a set 
of functions {x^} and has large values at energies where x^ (or {x^}) has large contributions 
in eigenstates (bands, molecular orbitals). PDOS can also be large in energy intervals with 
a large number of states. The weights in my PDOS are derived from Mulliken population 



analysis [42|, |43|. The ith eigenstate of the Kohn-Sham Hamiltonian ^jk is expanded in a 
finite basis 

V^ik(r) = ^C/,i(k)x^k(r), (4) 

where c^j(k) are expansion coefficients and x^tk are basis set functions - Bloch sums of 
equivalent atomic orbitals. The gross population of X/^k for the eigenstate ipi^ is 

GP.k(x,k) = IY1 {c,^ik)cUk)S,M + i(k)c..(k)5.,(k)) (5) 

V 

and PDOS for function x^ is 

PDOS^^{E) = J2Y1 GP,^ix,i.)LiE - E,k), (6) 

i k 

where Ei^ is Kohn-Sham eigenvalue corresponding to eigenstate ipik and L is a Lorentzian 
broadening function. 

OPWDOS is defined for two functions Xfi and Xv or between two sets of functions {{Xfi} 
and {Xi^}) and has large positive or negative values depending on whether the interaction 
between X/i and Xu ({X/i} and {Xu}) is bonding or anti-bonding, respectively. The use 



of these plots is demonstrated in Ref. |4J]. The overlap population of two orbitals and 



OPWDOS are defined as 






(7) 



III. RESULTS 



A. InP, GaAs, and InAs: Structural Parameters 

Fig. [1] shows energy level diagram for "spherically symmetric" In, As, Ga, and P atoms. 
Fig. [1] shows that group III A atomic species In and Ga have smaller ns — np and ns — {n — l)d 
energy spacings as compared to group VA species (As and P). The s — p energy spacings 
are 5.7-6.4 eV for In and Ga and 9.3-8.4 eV for As and P, respectively. 



InP, GaAs, and InAs are known to exist in several phases (See Refs. 10|, |45| and the 
references therein). The low pressure phase of all three (GaAs, InP, InAs) semiconductors is 
zinc-blende. While there had been some discussions on the nature of the first high pressure 
phase of GaAs, it is now generally agreed that the first transition takes place at approxi- 
mately 17 GPa and involves transition from zinc-blende (GaAs-I) to orthorhombic Cmcm 
phase (GaAs-II) 46] . The pressure release results in a transition to a cinnabar phase fol- 
lowed by a transition to the original zinc-blende phase 47|. No direct zinc-blende to cinnabar 
transition was observed. The first high-pressure phase in common "cation" InP and InAs 
was experimentally found to be NaCl phase 48N50| . The latter finding is supported by 



Density Functional calculations 51 



52|. 



In connection with InAs quantum dots in GaAs or InP matrix, the low pressure zinc- 
blende phase (Strukturbericht designation B3) is of primary interest. The structural param- 
eters of InP, GaAs, and InAs obtained from my Density Functional calculations with ZORA 
Relativistic Hamiltonian and APW+lo method are summarized in Table [B Table [T] also 
shows experimental values 45| which were measured at room temperature and the results of 
APW+lo calculations. Figure [2] shows the deviations between LSDA lattice constants and 
the experimental lattice constants. 

The finite temperature effects will increase the lattice constant. Once these temperature 
effects are taken into account, the "experimental" lattice constant is effectively reduced which 
will infiuence the conclusions about the accuracy of a given exchange-correlation functional. 



For example, in the case of GaAs, the temperature effects lead to the increase in the lattice 



constant by 0.3% from 5.638 A to 5.653 A 53|. Since LSDA underestimates 53|, |5J] the 



bond lengths, the inclusion of finite temperature effects into consideration will improve the 
agreement between the theory and experiment. 

Table [I] and Figure [2] show that LSDA ZORA relativistic and scalar relativistic results 
underestimate the lattice constants by, approximately, 0.6% (InP), 0.8 % (GaAs), and 0.5 
% (InAs). The consideration of finite temperature effects will further improve agreement 
between the theory and experiment and it is likely that the error of LSDA relativistic ZORA 
calculation for lattice constants is within 0.5%. 



In agreement with the established trend |55J, the relativity contracts the bond length. The 



DFT lattice constant decreases as the treatment of relativity changes from "nonrelativistic" 
(NREL) to fully relativistic (FREL). Variational treatment of spin-orbital coupling does not 
seem to affect the structural properties by much, the reduction in the equilibrium lattice 
constant upon going from SREL to FREL is very small (within 0.05%). It is important 
to include some kind of description of relativity for In - NREL lattice constants for InP 
and InAs are larger than experimental ones which contradicts to the established LSDA 



trend 



53 



Mi- 



Figure [3] shows LSDA bulk modulus calculated at different "levels of relativity" (FREL 
- fully relativistic with variational treatment of spin orbital coupling, SREL - scalar rela- 
tivistic, and NREL - nonrelativistic). My LSDA calculations reproduce the "stiffness" trend 
BinAs < Bjnp < BcaAs- The bulk modulus decreases upon going from NREL to SREL 
description and, then, slightly increases by going from SREL to FREL. In the case of InP 
and GaAs, the LSDA bulk modulus B^sda is smaller (less stiffer) than the experimental 
bulk modulus. This seems to contradict to the "established" LSDA trend of Blsda being 



too stiff, note, however, that by examining Table V of Ref. (53| one can not conclude this 
with respect to group IVA and group IIIA-VA semiconductors. In the case of InAs, ZORA 
FREL and SREL bulk modulus is by 2% more stiffer than the experimental bulk modulus. 
I find that my APW+lo calculations are in good agreement with the ZORA calculations. 
The APW+lo lattice constants are within 0.2% of the ZORA lattice constants. The bulk 
moduli obtained from the APW+lo calculations are smaller than those obtained with BAND 
program (FREL) by 2 %, 6%, and 15% for InP, GaAs, and InAs, respectively, but the 
"stiffness" trend BjnAs < BjnP < BcaAs is reproduced. 



B. InP, GaAs, and InAs: PDOS and OPWDOS Analysis 

Figs. mini and|8]show relativistic PDOS for InP, GaAs, and InAs, respectively. The PDOS 
is calculated at three values of the lattice constant corresponding to the tensile hydrostatic 
strain (Figs, (a) and (b)), equilibrium (Figs, (c) and (d)), and compressive hydrostatic strain 
(Figs, (e) and (f)). 

Let us consider PDOS at the equilibrium (experimental) lattice constants in InP, GaAs, 
and InAs (middle rows of Figs. HI |6l and [8]). The left (c) and right (d) columns show PDOS 
on the cation (In, Ga) and anion (As, P) atomic orbitals (AOs), respectively. The names 
"cation" and "anion" reflect the move away from covalent bonding and towards ionicity in 



InP, GaAs, and InAs semiconductors. The Hirshfeld charge analysis 56| performed in this 
work reveals that electron density transfers from regions near In and Ga ions into regions 
near As and P ions making In and Ga "positively" and As and P "negatively" charged, 
respectively. 

We find that PDOS of InP, GaAs, and InAs in the valence band energy region (up to 
16 eV below the Fermi energy, the Fermi energy is at zero) consists of four main "spectral 
features" . 

The first "spectral feature" in the PDOS is a peak of broad character just below the 
Fermi energy. The width of this peak is, approximately, 3.2 eV for InAs and InP and 4.4 
eV for GaAs. The main contribution to this "spectral feature" stems from anion and, to 
a lesser extent, from cation atomic orbitals (AOs) of p type. There is a small contribution 
from s and d cation orbitals and, in the case of GaAs and InAs, s and d anion orbitals. 

The second feature is a very sharp peak centered around 5.4 eV (InAs, InP) and 6.5 eV 
(GaAs) below the Fermi energy. The main contributions to this peak are cation AOs of s 
type. The smaller contributions to this peak are arising from anion and cation AOs of p 
type. The second feature has a slowly decaying tail which contributes to the very top of the 
valence band. 

The third peak in the PDOS is located, approximately, at 9.5 eV (InP), 12.2 eV (GaAs), 
and 10.9 eV (InAs) below the Fermi energy. The main contribution to this peak stem from 
anion s type orbitals as well as from cation p, s and d AOs. The last "spectral feature" is a 
very sharp peak in the PDOS due to the cation d orbitals and anion s type orbitals. 

The PDOS in the conduction band energy region is shown up to 9.5 eV above the Fermi 



8 



energy and consists of several closely spaced peaks. The main contributions in this energy 
window are the anion p and d and cation p, s, and d AOs. The PDOS does not show 
significant anion s-type contributions in the conduction band region. 

Figs. (a),(b) and (c),(d) demonstrate the effect of the hydrostatic tensile and compressive 
strains on the PDOS. The compressive hydrostatic strain 1) broadens PDOS peaks; 2) 
decreases magnitude of the PDOS peaks; 3) changes the relative position of the peaks. The 
"general structure" of the PDOS is preserved. 

Figs. El d and E] show OPWDOS for InP, GaAs, and InAs, respectively. The OPW- 
DOS is calculated at the experimental equilibrium lattice constants and provides a clear 
pictorial representation of orbital interactions in IIIA-VA InP, GaAs, and InAs zinc-blende 
semiconductors. 

The OPWDOS in the valence band energy region is of mostly bonding character. The 
OPWDOS shows that the very top of the valence band is strongly stabilized by cation-anion 
p — p orbital interactions. The energy region corresponding to the second PDOS peak is 
characterized by the cation s - anion p and cation p - anion s bonding interactions as well 
as a slightly anti-bonding cation s - anion s interactions. The energy region corresponding 
to the third PDOS peak is characterized by the cation p - anion s interaction. 

The OPWDOS in the conduction band energy region is mostly of the anti-bonding char- 
acter and has peaks that are generally higher in magnitude than the OPWDOS peaks in 
the valence band region. The latter is a demonstration of a well known fact that the "anti- 
bonding is more anti-bonding than bonding is bonding" (see Ref. [57| and the references 
therein). The bottom of the conduction band is strongly "destabilized" by cation-anion 
s — p, p — s, and s — s interactions. It is interesting to note that the cation-anion p — p 
interaction for the very bottom of the conduction band (within 1 eV) is bonding. At the 
higher energies, the p — p interaction becomes strongly anti-bonding. 

C. InP, GaAs, and InAs: Electronic Structure Parameters 

The PDOS and OPWDOS presented in Section IIIIBI describe the electronic structure 
qualitatively. Tables [IT] and IIVI show results for the parameters of the electronic structure 
for InP, GaAs, and InAs IIIA-VA zinc-blende semiconductors. These parameters are the 
band gaps for the transitions between the high-symmetry points of the Brillouin zone E^^^ = 



E2 — El, split-off energies at the F and L points of the Brillouin zone, position (with respect 
to the Fermi energy) and width of the cation d-band at the F point, and the width of 
the upper part of the valence band. Table IIVI summarizes the relative volume deformation 
potentials for a specific 1—7-2 transition defined as 

1^2 _ dE2i , . 

The electronic band structure parameters are computed at several levels of relativity (FREL 
- fully relativistic, SREL - scalar relativistic, and NREL - nonrelativistic) as well as with 
the APW+lo method (fully relativistic approach). The results of these calculations are 
compared with available experimental data. 

In agreement with previous work. Table [TTl shows that LSDA KS DFT severely underes- 
timates band gaps. The small LSDA band gap becomes even smaller when one introduces 
explicit description of relativity. For example, in the case of GaAs, the direct gap -E'J^"^^'^ 
decreases from 1.0 eV to 0.4 eV as the "level of relativity" changes from nonrelativistic to 
fully relativistic. The relativity especially strongly affects band gaps at F and L points of 
the Brillouin zone. 

In the Table llllt I summarized orbital populations for three top valence bands (split-off 
Tjy, light-hole, and heavy hole Fg?,) and the lowest conduction band (Fgc) at the high- 
symmetry points of the Brillouin zone. From Table IIIH one can see that the bottom of 
the conduction band at the F and L points has significant contributions from cation s-type 
AOs, whereas at the X point the bottom of the conduction band is made up of cation p-type 
orbitals. The reason for the strong "relativistic" band gap decrease at the F and L points is 
that the conduction s states are stabilized by relativity stronger than the valence p states. 



The stabilization itself stems from the relativistic "contraction" of the atomic orbitals 58|. 
Overall, the gaps in the relativistic description decrease substantially which might affect the 
conclusions with respect to the accuracy of a given exchange-correlational functional for the 
band gap calculation. 

For a fixed lattice constant, LSDA also poorly describes the energy differences within the 
conduction band. For example, E{Xc) — E{Lc) energy difference is 508 meV from ZORA 
FREL calculations, whereas the experimental value is 160 meV. The relative position of the 
conduction band minima is described only qualitatively E^^^^'^ < E^^^^^'^ < E^^^^^ . 

The strong underestimation of the band gaps may result in the wrong energetic order of 
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bands in some specific points of tlie Brillouin zone. Tliis liappens at tlie F point for InAs at 
the equilibrium lattice constant, where a conduction band Fqc is strongly stabilized and lies 
below the split-off band Tj^. The strain also may affect the energetic order of bands. 

Figures [10] and [TT] show scalar relativistic and fully relativistic band structure plots for 
InP, GaAs, and InAs, respectively. The band structure is computed along the edges con- 
necting the high-symmetry points of the Brillouin zone. The Cartesian coordinates of these 
high-symmetry points are summarized in Table |Vl From these band structure plots one can 
see the band crossing at the F-point for InAs at the equilibrium lattice constant. This band 
crossing also occurs in InP and GaAs subjected to the tensile hydrostatic strain. 

It is found that the split-off energies at the F and L points are reproduced very accurately 
within ZORA fully relativistic approach, the agreement with the experiment is a few % or 
several me Vs. The fully relativistic treatment for the calculation of the split-off energies is 
essential as both scalar relativistic/nonrelativistic calculations lead to the six-fold (including 
spin) degeneracy of the valence band at the F-point (Fig. [TOl) . 

The relative volume deformation potential describes how fast a given band gap changes 
with volume. The negative (positive) relative deformation potential (in our definition) means 
that the band gap increases (decreases) as volume decreases. Table [TTl shows that for both 
Tv — )■ Fc and Fy — )■ Lc transitions, both gaps are increasing as volume decreases (negative 
deformation potential), whereas for the Ty — > Xq transition, the band gap decreases (the 
deformation potential is positive). The differences in the sign of the deformation potential 
for Ty — !■ Fc and Fy — )■ Lc transitions on one hand and Fy — )■ Xc transition on the other 
hand are attributed to the "different nature" of the conduction band minimum at these 
points (see Table UTTl) . 

The calculated absolute magnitude of the "rate of change" in the gap is the largest for 
GaAs and decreases for semiconductors with a larger lattice constant (InP, InAs). The ex- 
perimental relative deformation potentials are obtained from the direct band gap pressure 
dependence coefficients and experimental bulk moduli. The experimental trend in the ab- 
solute magnitude of the "rate of change" in the gap is GaAs, InAs, and InP. Note, however, 
that experimental uncertainties for the relative deformation potential can be as large as 1 
eV. The relative deformation potentials are quite close for fully relativistic and scalar rela- 
tivistic calculations and, therefore, a fully relativistic calculation of this quantity does not 
seem essential, at least, for the transitions considered in this work. 
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Finally, there are two other quantities which sensitively depend on fully relativistic calcu- 
lation - the cation rf-band width and the width of the upper part of the valence band UVBW. 
Both band widths increase as the treatment of relativity changes from non-relativistic to 
fully relativistic level. For UVBW, the increase is 6%-8% (0.4-0.5 eV). The rf-band width 
increases dramatically from 0.10-0.15 eV to 0.9 eV. 

The agreement between ZORA fully relativistic and APW-l-lo calculations is exceptionally 
good, especially, for the relative deformation potentials, E^, SE^, and for UVBW. The split- 
off energies are reproduced within 10 meVs, and the gaps usually agree within 20 me Vs. 

IV. CONCLUSIONS 

First-principles full potential calculations based on the Zeroth Order Regular Approxi- 
mation (ZORA) relativistic Hamiltonian and the Kohn-Sham form of Density Functional 
Theory (KS DFT) were reported for group IIIA-VA (InAs, GaAs, InP) semiconductors. 
The effects of relativity were elucidated by performing fully relativistic, scalar relativistic, 
and nonrelativistic calculations. The inclusion of relativity led to the contraction of the 
bond length, strong stabilization of the conduction band at F and L points of the Brillouin 
zone, and broadening of the upper part of the valence band. The inclusion of relativity at 
least on the scalar relativistic level was found to be essential for the accurate calculation 
of structural properties. Electronic band structure parameters were determined including 
split-off energies, band gaps, deformation potentials, and populations at the high-symmetry 
points of the Brillouin zone. It was found that the split-off energies can be determined 
with very good accuracy. In agreement with previous work, LSDA KS DFT severely un- 
derestimates band gaps which may result in the wrong energetic order of bands at specific 
points of the reciprocal space. It was found that relativistic LSDA describes the sequence 
of conduction band minima only qualitatively. The relative band gap deformation poten- 
tials were determined and compared with the available experimental data. The relativistic 
LSDA relative deformation potentials at the F point were found to be too small. The nature 
of the chemical bonding at the equilibrium and under hydrostatic strain was investigated 
using projected (PDOS) and overlap population weighted density of states (OPWDOS). It 
was found that OPWDOS in the valence and conduction band energy regions is of mostly 
bonding and anti-bonding type, respectively. ZORA results were compared with Augmented 
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Plane Wave plus Local Orbitals method (APW+lo) and a good agreement between the two 
sets of calculations was established . Viability and robustness of the ZORA Hamiltonian for 
the investigation of electronic and structural properties of semiconductors was reaffirmed. 
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FIG. 1: "Spherical atom" energy levels. 
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FIG. 2: LSDA error (%) between KS DFT calculated lattice constant and experimental lattice 
constants at different "levels of relativity" (FREL - fully relativistic, SREL - scalar relativistic, 
and NREL - nonrelativistic). LSDA lattice constants (FREL, SREL) underestimate experimental 
constants. In the case of InAs and InP, the NREL DFT lattice constant is larger than the exper- 
imental one which contradicts to the established trend. The experimental lattice constants were 
determined at the room temperature. 
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FIG. 3: LSDA bulk modulus calculated at different "levels of relativity" (FREL - fully relativis- 
tic, SREL - scalar relativistic, and NREL - nonrelativistic) . The LSDA calculations (this work) 
reproduce the "stiffness" trend BjnAs < BjnP < BcaAs- 
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4: PDOS in InP. Results of fully relativistic calculations for several lattice constants. 
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FIG. 5: OPWDOS in InP (a = 5.869 A). Results of scalar relativistic calculations. 
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PDOS in GaAs. Results of fully relativistic calculations for several lattice constants. 
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FIG. 7: OPWDOS in GaAs (a = 5.648 A). Results of scalar relativistic calculations. 
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FIG. 8: PDOS in InAs. Results of fully relativistic calculations for several lattice constants. 



c 



O 

Q 
a 



0.14 

0.12 
0.10 
0.08 
0.06 

0.04 
0.02 
0.00 



3 

(A 
O 

a 

a. 



c 
3 



to 
O 
O 

CL 



-^ — T — I — ■r^ — I — 1 — K — » — t — t — 1 — T — r — 1 — T — * — r-r 



Ind 


- 




Inp ^ 


^H : 




ins H 


^^m ~ 






' 


tn 




_ 


■e 




" 


c 




- 


3 




< 


.d 




: 


* 


1 


^ 


(O 


.1 


h 


O 


1 11 


. 


Q 


1 u 


i 


0- 



I I ^^> I ■ ■ ■ 




-15.0 -10.0 -5,0 0.0 
energy [eV] 

(a) fl=6.400 A 



5.0 



0.14 
0.12 

0.10 
0.08 
0.06 
0.04 
0,02 
0,00 



Ind 
In p 
Ins 



•15.0 -10.0 -5.0 0,0 
energy [eV] 

(c) a=6.058 A 



5,0 



0,14 
0,12 
0,10 
0.08 
0.06 
0.04 
0.02 
0.00 



-15.0 -100 -5.0 0.0 
energy [eV] 

(e) «=5.600 A 




Ind 






inp ^ 


^B = 




Ins ^ 


^^ - 








in 






ts 






C 






3 






£i 






to 






CO 






o 






Q 






Q. 




3 

o 

Q 
a. 




-15.0 -10,0 -5.0 



energy [eV] 

(b) a^6.400 A 




-15.0 -10,0 -5,0 0,0 5.0 



energy [bV] 

(d) a=6.058 A 



0,14 
0,12 
0,10 
0,08 
0,06 
0.04 
0,02 
0.00 



[ ' ' ' ' I 
As d 
As p ^ 
Ass m 




-15.0 -10,0 -5.0 0-0 
energy [eV] 

(f) fl=5.600 A 



50 



23 



FIG. 9: OPWDOS in InAs (a = 6.058 A). Results of scalar relativistic calculations. 
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FIG. 10: Band Structure of InP, GaAs, and InAs. Results of scalar relativistic calculations. 
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FIG. 11: Band Structure of InP, GaAs, and InAs. Results of fully relativistic calculations. 
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TABLE I: 

Structural parameters of InP, GaAs and InAs IIIA-VA semiconductors (zinc-blende crystal struc- 
ture). LSDA equilibrium lattice constant a^g (A), bulk moduli B (kbar=0.1GPa) and pressure 
derivatives B'. Experimental values are taken from Ref. [43]. The experimental bulk modulus is 
computed from the elastic constants cn and ci2 as B = (cn + 2ci2)/3. 

Method Ueq (A) B (kbar) B' 

InP 

ZORAFREL 5.831 699 4.5 

ZORA SREL 5.833 697 4.5 

NREL 5.873 709 4.4 

APW+lo 5.823 685 4.4 



EXP. 5.869'' 710' 



4 



GaAs 



ZORAFREL 5.607 736 4.4 

ZORA SREL 5.610 732 4.4 

NREL 5.627 745 4.5 

APW+lo 5.610 691 4.4 

b 7Ka<i 



EXP. 5.653" 756' 



InAs 



ZORAFREL 6.026 590 4.7 

ZORA SREL 6.029 590 4.8 

NREL 6.069 614 4.4 

APW+lo 6.033 511 4.6 



EXP. 6.058^= 579' 



4 



Experimental temperature: "T=291.15 K, T=298.15 K, '^T=300 K, room temperature 
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TABLE II: 

Parameters of electronic structure for InP, GaAs, and InAs IIIA-VA semiconductors (zinc-blende 
crystal structure). All energies are in eV. Egap - energy gaps for transition between high-symmetry 
points of the Brillouin zone, A^q and A^q - magnitude of spin-orbital splitting at the top of the 
valence band at T and L points, E^ - d band position at the F point, 5E^ - d band width at the 
r point, and UVBW - the upper valence band width. 



Method 


gap 


gap 


gap 

InP 


Ago 




E', 


^El 


UVBW 


ZORA FREL 


0.437 


1.606 


1.292 


0.108 


0.120 


-14.19 


0.89 


5.89 


ZORA SREL 


0.475 


1.644 


1.330 


0.0 


0.0 


-14.09 


0.16 


5.85 


NREL 


1.014 


1.590 


1.781 


0.0 


0.0 


-15.07 


0.14 


5.51 


APW+lo 


0.416 


1.575 


1.257 


0.100 


0.112 


-14.21 


0.89 


5.90 


EXP. 


1.420 


2.384 


2.014 


0.110 


0.130 














GaAs 












ZORA FREL 


0.196 


1.257 


0.748 


0.349 


0.213 


-15.01 


0.48 


7.00 


ZORA SREL 


0.313 


1.372 


0.865 


0.0 


0.0 


-14.85 


0.08 


6.89 


NREL 


0.965 


1.418 


1.197 


0.0 


0.0 


-15.38 


0.08 


6.69 


APW+lo 


0.173 


1.240 


0.735 


0.341 


0.209 


-15.02 


0.48 


6.99 


EXP. 


1.520 


1.979 


1.819 


0.340 


0.220 














InAs 












ZORA FREL 


-0.577 


1.325 


0.661 


0.363 


0.262 


-14.46 


0.88 


6.02 


ZORA SREL 


-0.459 


1.444 


0.784 


0.0 


0.0 


-14.27 


0.13 


5.90 


NREL 


0.405 


1.453 


1.337 


0.0 


0.0 


-15.26 


0.11 


5.49 


APW+lo 


-0.594 


1.325 


0.653 


0.351 


0.250 


-14.48 


0.88 


6.03 


EXP. 


0.420 


2.24 


1.133 


0.380 


0.267 
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TABLE III: 

Orbital populations of the band edges in InP, GaAs, and InAs semiconductors (zinc-blende crystal 
structure) in high-symmetry points of the Brillouin zone at the equilibrium lattice constants. SO - 

split-off band, LH - light hole band, HH - heavy hole band, and CB - conduction band minimum. 

ZORA FREL LSDA calculation. 

InP 
T X L 

In.s In.p P.s P.p In.s In.p P.s P.p In.s In.p P.s P.p 

SO 0.000 0.112 0.000 0.809 0.619 0.000 0.000 0.327 0.451 0.328 0.042 0.179 

LH 0.000 0.072 0.000 0.878 0.000 0.360 0.000 0.630 0.000 0.237 0.000 0.752 

HH 0.000 0.072 0.000 0.878 0.000 0.339 0.000 0.645 0.000 0.209 0.000 0.786 

CB 0.920 0.000 0.067 0.000 0.000 0.621 0.030 0.000 0.430 0.252 0.014 0.234 

GaAs 

r X L 

Ga.s Ga.p As.s As.p Ga.s Ga.p As.s As.p Ga.s Ga.p As.s As.p 

SO 0.000 0.118 0.000 0.755 0.523 0.000 0.000 0.402 0.452 0.352 0.015 0.183 

LH 0.000 0.128 0.000 0.738 0.000 0.435 0.000 0.539 0.000 0.294 0.000 0.661 

HH 0.000 0.128 0.000 0.738 0.000 0.453 0.000 0.524 0.000 0.279 0.000 0.636 

CB 0.728 0.000 0.216 0.000 0.000 0.556 0.040 0.000 0.305 0.298 0.092 0.130 

InAs 
T X L 

In.s In.p As.s As.p In.s In.p As.s As.p In.s In.p As.s As.p 

SO 0.000 0.100 0.000 0.829 0.610 0.000 0.000 0.303 0.485 0.279 0.019 0.217 

LH 0.000 0.081 0.000 0.843 0.000 0.339 0.000 0.637 0.000 0.242 0.000 0.750 

HH 0.000 0.081 0.000 0.843 0.000 0.360 0.000 0.611 0.000 0.210 0.000 0.780 

CB 0.897 0.000 0.093 0.000 0.000 0.605 0.039 0.000 0.354 0.347 0.026 0.143 
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TABLE IV: 

Relative volume deformation potentials ay (eV) for InP, GaAs, and InAs IIIA-VA semiconductors 
(zinc-blende crystal structure) for specific transitions. The different signs of the deformation poten- 
tial is attributed to the "different nature" of the conduction band minimum. Experimental volume 
deformation potentials are obtained from the direct band gap pressure dependence coefficient and 
bulk modulus. 



Method fly''^^'' alX^^'' alX^^'' 



InP 



ZORAFREL -5.44 1.64 -2.42 

ZORA SREL -5.45 1.64 -2.42 

NREL -5.12 1.60 -2.65 

APW+lo -5.44 1.64 -2.40 



EXP. -5.7 



GaAs 

ZORAFREL -7.46 1.83 -2.78 

ZORA SREL -7.52 1.79 -2.82 

NREL -7.79 1.81 -3.06 

APW+lo -7.41 1.84 -2.74 

EXP. -8.0,-9.2 

InAs 

ZORAFREL -5.04 1.59 -2.11 

ZORA SREL -5.09 1.56 -2.15 

NREL -5.11 1.51 -2.58 

APW+lo -5.01 1.59 -2.09 

EXP. -6.6 
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TABLE V: 

Cartesian coordinates of special k-points used to construct the band structure plots (special k- 
points are shown with vertical lines), a - lattice constant. 

Symbol Coordinates 

^[1,-1/4,1/4] 
L f [1/2,-1/2,1/2] 

r f [0,0,0] 

^[1,-1/4,1/4] 
X f [1,0,0] 

r f [0,0,0] 

W ^[1,0,1/2] 

^[1,-1/4,1/4] 

r vfO'O'O] 
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